Nonchaotic Stagnant Motion in a Marginal Quasiperiodic Gradient System 
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A one- dimensional dynamical system with a marginal quasiperiodic gradient is presented as a 
mathematical extension of a nonuniform oscillator. The system exhibits a nonchaotic stagnant mo- 
tion, which is reminiscent of intermittent chaos. In fact, the density function of residence times 
near stagnation points obeys an inverse-square law, due to a mechanism similar to type-I intermit- 
tency. However, unlike intermittent chaos, in which the alternation between long stagnant phases 
and rapid moving phases occurs in a random manner, here the alternation occurs in a quasiperiodic 
manner. In particular, in case of a gradient with the golden ratio, the renewal of the largest res- 
idence time occurs at positions corresponding to the Fibonacci sequence. Finally, the asymptotic 
long-time behavior, in the form of a nested logarithm, is theoretically derived. Compared with the 
Pomeau-Manneville intermittency, a significant difference in the relaxation property of the long-time 
average of the dynamical variable is found. 

PACS numbers: 05.45.Ac, 05.45.Pq 



I. INTRODUCTION 

Temporal intermittency is the occurrence of a signal 
accompanied by random alternation between long lam- 
inar phases and relatively short bursts, which is widely 
observed in nonequilibrium dynamical systems. For ex- 
ample, the Rayleigh-Benard convection [ij, the Belousov- 
Zhabotinsky reaction f3|, and an rf-driven Josephson 
junction [3|] exhibit intermittent phenomena. In the field 
of nonlinear physics, intermittency is tacitly understood 
as the occurrence of intermittent chaos. 

Here we confine ourselves to the Pomeau-Manneville 
intermittency Q . The onset of intermittent chaos is as- 
sociated with a loss of stability of periodic motion, which 
is classified into three types: type-I (saddle-node bifur- 
cation) , type-II (Hopf bifurcation) , and type-Ill (subhar- 
monic bifurcations). The Pomeau-Manneville (PM) sys- 
tem, Xn+i — T{xn) = Xn + ax^ (mod 1) (z > 1, a > 0), is 
a typical model of a Poincare map for intermittent chaos 
The system is called non-hyperbolic, since it has 
an indifferent fixed point a; = with T'{x) = 1, in the 
neighborhood of which laminar motions are generated, 
and the ergodic measure pix) localizes as p(x) oc 

One of the important problems in such a non- 
hyperbolic system is the appearance of nonstationarity. 
In a nonstationary regime where z > 2, the ergodic mea- 
sure p{x) cannot be normalized, and the f~'^ {v > 1) 
fluctuation is generated. Nonstationarity can be shown, 
for example, in the behavior of the renewal function 
H{n) = E{Nn), which is the ensemble average of the 
number Nn of chaotic bursts during the time interval 
(0, n] [1, H, Therefore, the renewal rate defined by 
H{n)/n represents the average occurrence probability of 
the bursts during (0, n]. In stationary regimes, this prob- 
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ability does not depend on time, H(n)/n ~ l/r. How- 
ever, in the nonstationary regime (z > 2), it depends on 

_ 2 — 2 

time as H{n)/n oc n . In particular, at the critical 
point (z = 2), it behaves as H{n)/n oc 1/lnn. Thus, 
the renewal rate indicates the nonstationarity of the in- 
termittent chaos. 

Although many studies of intermittency have dealt 
with chaotic systems, intermittent dynamics are also 
observed in nonchaotic systems. Intermittent strange 
nonchaotic attractors (SNAs) are generically created in 
quasiperiodically forced systems through quasiperiodic 
saddle-node bifurcations [lOl [TlJliL through quasiperi- 
odic subharmonic bifurcations ij l3] , and through sev- 
eral types of crisis [llf, where the scaling behavior is 
characteristic of type-I, type-Ill, and crisis-induced inter- 
mittency, respectively. Unlike the analogy between the 
intermittencies in SNAs and those in chaotic systems, to 
the best of our knowledge, the difference is not clear, ex- 
cept for the original difference in the sign of the largest 
nontrivial Lyapunov exponent [l^, El, [13, [H, [3| • 

In this study, we present a marginal quasiperiodic gra- 
dient system (MQPGS) as a mathematical extension of a 
nonuniform oscillator x = 1 — Acosx, which arises in sev- 
eral fields, such as electronicSjCondensed-matter physics, 
mechanics, and biology [H, [l^, and study the stagnant 
motion generated by the system. The term "nonchaotic 
stagnant motion" in this paper represents anomalous dy- 
namics accompanied by long laminar phases and those 
interruptions not based on chaotic dynamics. Stagnant 
motion in a MQPGS is related to the spatial quasiperi- 
odicity of the gradient, which differs from the temporal 
quasiperiodicity in intermittent SNAs. The quasiperiodic 
gradient could be implemented in some Josephson junc- 
tion circuits. Our main result is the theoretical deriva- 
tion of the asymptotic divergence of the displacement in 
a MQPGS. By comparing the asymptotic behavior in the 
MQPGS with the renewal rate in a PM system, we clar- 
ify a significant difference in the relaxation property of 
the long-time average of the dynamical variable. 
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This paper is organized as follows: Section II intro- 
duces the MQPGS and mentions some implementation 
methods, Section III presents the analysis of a two-period 
case of a MQPGS, and Section IV gives the derivation of 
the density function of residence times near stagnation 
points. Section V clarifies the parameter dependence of 
the largest residence time from a number-theoretic point 
of view. Section VI gives the derivation of the asymptotic 
behavior of the MQPGS, and the last section presents a 
comparison between the nonchaotic stagnant motion of 
MQPGS and the intermittent chaos of the PM system 
from the viewpoint of asymptotic behavior. 



II. MARGINAL QUASIPERIODIC GRADIENT 
SYSTEM 



A nonuniform oscillator 
tion 



15| is described by the equa- 



X 



1 — A cos X 



(1) 



in the time scale normalizing the phase-averaged angular 
frequency. Here A is the control parameter, and A = Q 
corresponds to the uniform oscillator. The nonuniform 
oscillator has a phase-drift (A < 1) or phase-locked state 
[A > 1) as a result of saddle-node bifurcation. There 
are many oscillatory phenomena explained by Eq. ([1]), 
such as those in oscillating neurons, firefly flashing, a 
Josephson junction and an overdamped pendulum driven 
by a constant torque [l^, [lB| ■ 

As a mathematical extension, it is natural to ask what 
happens when the right-hand side of Eq. ([T]) becomes 
quasiperiodic. Thus, we introduce an MQPGS described 
by a one-dimensional ordinary differential equation 



dx(t) 
dt 



\ — Ai cos{kix + 5i) — A2 cos{k2X + 62), (2) 



where Ai > 0, k^, Si (i = 1,2) are parameters, fci/fc2 is 
irrational, and Ai + A2 = 1. We are also interested in 
the rational system (fci/fc2 is rational), since the MQPGS 
is indistinguishable from well-approximated rational sys- 
tems in a finite time and with finite resolution. 

The MQPGS could be implemented in asymmetric 
multijunction superconducting quantum interference de- 
vices (SQUID) modeled after 3JJ SQUID ratchet pro- 
posed by Zapata et al. [l3|. Figure 1 shows the schematic 
representation of an asymmetric multijunction SQUID 
threaded by a flux <I>oxt , where the left and right branches 
contain m and n identical Josephson junctions, respec- 
tively. The SQUID is driven by a dc current Jdc, which 
splits into two branch currents, Ji (left) and Jr (right). 
We assume that the junctions are described b y th e re- 
sistively shunted junction (RSJ) model [H, IT9ll20| with 
resistance R and critical current /. Furthermore, we con- 
sider the overdamped limit, {2e/h)IB?C -C 1, where the 
junction capacitance C is negligible (see Fig. 1(b)), and 
also the Nyquist noise current is neglected. Then, the 
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FIG. 1: (a) Asymmetric multijunction SQUID for an imple- 
mentation of MQPGS. (b) Resistively shunted junction model 
of a single Josephson junction in overdamped limit. 



phase 4'ii, across the single junction i on the left branch, 
obeys the following equation: 



2eR 



(jjii + Ii sin((/);i) = Ji 



(3) 



On the right arm, the phase (j>rj of the junction j obeys 
the equation obtained by replacing labels I and i in 
Eq. ^ with r and j, respectively. If identical initial 
conditions are assumed, the only solutions for each junc- 
tion are phase-locked states, (jjn = 4>i2 = ■ ■ ■ = (pi^ = (pi 
and (p^-i — (f>r2 = ■ • • = 4'rn = 4>r- the limit where the 
total loop inductance L is negligible, \LI\ <^ $0i the to- 
tal flux is approximately equal to the external flux $cxt, 
where $0 = ft/2e is the flux quantum. Then, the integra- 
tion of the gauge invariant phase around the loop yields 
mpi — ncpr — — 27r<I>oxt/*I'o + Stts (s € N). Consequently, 
the Kirchhoff's current law Jdc = Ji + Jr is given by 



TO 

2eR K^^l^ 



Jdc - Ii sin( 



X sm I — (f>i + 277 

n n 



which approximately describes the MQPGS with ^2 /fci ~ 
m/n, under the condition of J^c = Ii + Ir: if all the 
variables are transformed properly. Note that the better 
approximation of fc2/fci ~ m/n yields the better imple- 
mentation. 

Also, the MQPGS could be implemented in a dc- 
driven circuit consisting of overdamped Josephson junc- 
tions coupled by ideal transformers. However, we omit a 
detailed description, since this is not the purpose of this 
paper. 



III. NONCHAOTIC STAGNANT MOTION 

In the following sections, we analyze dynamics of 
MQPGS, focusing on a typical case: 



dx{t) 



1 cos(27ra;) cos(27rfca;), (4) 

dt 2 ^ ' 2 ^ ^ ^ ^ 



where fc is a control parameter, and x(0) 



Xq. 



The 
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TABLE I: Parameters for lines in Fig. 2 
Line k Continued fraction " Approximation 



(b) 

V"/ 610 



[0; 1", 2, p5] 
[0; P^ 10, 
[0; 1"] 



0.618033988749894 
0.618033617353155 
0.618033260658830 
0.618032786885245 



"For example, the notation [0; l'^*, 2, 1^^] means a finite con- 
tinued fraction [go; 91, 92, ...,94o] with partial quotients, go = 0, 
91 = 92 = ■ ■ ■ = gi4 = 1, 915 = 2, and qig = ■ ■ ■ = q^g = 1. 
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FIG. 2: Displacement trajectories x{t) starting from the same 
initial condition xq = 1/2 for various values of parameter k 
near the inverse golden ratio tp^^ = ■ k values for each 
line are given in Table 1. The inset shows the velocity v{t) 
for k — (fi~^ ■ 



displacement x{t) never decreases but frequently slows 
down due to the quasiperiodic gradient. Figure 2 shows 
the trajectories of the displacement x{t) starting from 
the same initial condition = 1/2 for various values of 
k near the inverse golden ratio ip^^ — ■ In the case 
where fc is a rational number, p/q, with co-prime integers 
p and q, the displacement x(t) finally goes to one of the 
equilibrium states x = qn{n € Z). Actually, the line (d) 
in Fig. 2 is for k = 377/610, which has an equilibrium 
point ai X — 610. If the parameter k is p' /q' , which dif- 
fers slightly from p/q, and q' > q, the displacement x{t) 
takes a long time to pass through the vicinities of equilib- 
rium points for k = p/q. We call these stagnant phases. 
Lines (a)-(c) in Fig. 2 represent the displacements for 
several fc-values slightly different from 337/610. We ob- 
serve the alternation between long stagnant phases and 
rapid moving phases, which is reminiscent of intermit- 
tent chaos (see inset of Fig. 2). When k is irrational, the 
equilibrium points vanish except for the origin, but the 



motion stagnates in the vicinity of a; = gn if the value of 
k is close to p/q. Note that every rational A: is a bifurca- 
tion point; i.e., the system is structurally unstable with 
respect to k. 

When the motion stagnates, the factor cos(27ra;) must 
approach unity. Thus, the coordinates of the stagnation 
points are restricted to the vicinities of integers n. Hence, 
we consider the relative dynamics of the variable y 
x — n (— ^ < y < 5) within each cells /„ = [n — ^, n 
and introduce new parameters: 

_ J A:n (mod 1) if A:n (mod 1) < 1/2, 
^" " \ fcn(mod 1) - 1 if fcn(mod 1) > 1/2, 

Sn = |7n| = min[fcn (mod 1), I — kn (mod 1)]. 
Then, the dynamics in the n-th cell is written as 

^ - 1 - i cos(27ry) - i cos(27r(7„ + ky)). (5) 

The parameters 7„ and e„ are uniquely determined for 
each n in — i < 7„ < i and in < e„ < i, respectively. 
They relates to the amount of the stagnation in each cell. 



IV. DENSITY FUNCTION OF RESIDENCE 
TIMES 

In the following, let k be irrational, and the initial 
condition be xq = 1/2. Then, the displacement x{t) 
visits every cell In {n > 1) only once in the course of 
time. Let us define the residence time Tn as the pe- 
riod for which the object stays in the n-th cell /„ as 

Tn = f^i -, — I jT. — r-^ — T^w — rr^- Due to the sym- 

metry of the integral interval, 7„ in the integral can be 
replaced with £„. Thus, the residence time T„ is written 
as a function of e„. 



T — 

-^n — 



dy 



cos(27ry) — i cos(27r(e„ + ky)) 



(6) 



To estimate the residence time, consider the narrow y- 
region for each cell, in which cos(7ry) and cos(7r(7„ + ky)) 
approach unity simultaneously, and are approximated by 
the second-order Taylor expansion around zero for each 
phase. When fc < 1 is satisfied, at most only one such 
region is present within each cell, and is expressed by 
parameter r near unity as follows: 

In,r = { 2/ 1 cos(27ry) > T A cos(27r(7„ -I- ky)) > r}, 

= (max [-5{r), -{6{r) +7„)/fc] ,mm[6{r), {6{r) - 7„)/fc]) 

where 6{r) — ^arccos(r). For fc > 1, Eq. ([5]) is reduced 
to a similar form for < fc' < 1 using scale transforma- 
tion fc' — 1/fc, x' — kx, and t' = kt. Therefore, we limit 
our discussion to the case of < fc < 1. In each region 
In,r, Eq. (|4]) is approximated by the second-order Taylor 
expansion. 



^^.2(l-ffc2)(y + 



l-|-fc2 



l-|-fc2' 



(7) 



4 



The second term of Eq. ([7]) gives the minimum velocity given by the integration of Eq. ([7|), 
in each cell. The residence time T„ in the region /„ ^ is 



J 



1 



7r2(l + fc2) 



(en^O). 

7r£„ 



■ arctan 



min[5(r),(5(r-)-7.)/fc] + -^a 



,[_5(r),-(5(r)+7„)/fc] + ^ 



(8) 



In the complementary regions /„ \ In,r, the residence 
times are shorter than 1/(1 — r), since y > 1 — r. There- 
fore, the residence time T„ is dominated by Tn^r for small 



(9) 



This is the same as the universal scaling of type-I inter- 



mittency with an index of —1/2, since r„ ~ (Tr^e^)"^/^ 
when the channel width is Tr^e'^/{1 + fc^). 

After sufficient displacement x{t), the density func- 
tion of observed e„ converges to a uniform density 
F{en) = 2 on the interval (O, ^), since the values of kn 
(mod 1) (n € N) are distributed uniformly on the torus 
interval (0, 1) [2]|. As a result, the density function P{T) 
of residence times is obtained by P{T) = —F{e)^ and 
shows an inverse-square law when the residence time T 
is sufficiently large, 

^(^)^^ (T^^). (10) 

Note that this density function does not have any finite 
moments. 

The numerical results agree with Eqs. ^ and pn]) in 
their asymptotic regions, as shown in Fig. 3. It is signifi- 
cant that the asymptotic behavior of the density function 
P{T) is independent of parameter fc, provided that k is 
irrational. Nevertheless, the observation time for obtain- 
ing the universal density function depends substantially 
on k. 



V. RENEWAL PROCESS OF THE LARGEST 
RESIDENCE TIME 

Let us consider that the displacement x{t) has passed 
through the finite space interval [xq, ni + i), where ni 
is the index of the latest cell /„, that the displacement 
passed through. Then, the largest residence time during 
the passage up to the cell /„, is defined by 



n<ni 



-'-rKni + 1 



max{T„ \1 < n < ni}. 



When the number n; of passed cells increases, the largest 
residence time T^^^^ is renewed. This section clarifies 
the positions of cells where the renewal occurs. 



P{T) 



0.01 



0.001 



0.0001 
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FIG. 3: The density function P{T) of residence times for 
k = tp~^ ■ Residence times from n = 1 to 5000 are used. The 
bin width is set to 1. The solid line represents a theoretical 
estimate obtained from Eq. (|10p. Deviations are observed for 
large T, where the statistics are poor. The inset shows the 
residence times T„ as a function of e„. The residence times 
Tn from n = 1 to 2000 are plotted for k = (fi'^ ■ The solid line 
is the theoretical estimate from Eq. ((O}. Deviations from the 

theoretical estimate \Tn — \ converge to about 0.063±0.001 

when e„ < 0.001. 



The parameter k can be expressed by the 
continued fraction representation as A; = 
bo; li, 92, • • • di] , where go G Z, q, e N, 0, e 

E, and Ot > 1. {qi} are called partial quotients, defined 
by the formula qi = [9i\ . 9i is the z-th complete quotient, 
generated by the recursion formula 0i — Qi + j^- Replac- 
ing 9i in the continued fraction with q^, we get the z-th 
convergent of k, i.e., Pt/Qi [go; 9i, 92, ••• ,9j-i, Qt]- 
These successive convergents are generated by the 
following recursion relations: 

P, = g,P,_i + P,_2, P-2 = 0, P-1 = 1, 

Q^ = q^Q^-l + Q.-2, Q -2 = 1, Q-1 = 0. (11) 

Note that the sequences {Pi} and {Qi} increase mono- 
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tonically. 

Let m be the nearest integer to kn. Then, £„ is given 
by Sn — \kn — m\. Using the above parameters, we can 
rewrite e„ for n — Qi (i > I) as 

eQ.^m.-P.\ = ^ {^>l), (12) 

which is derived in Appendix A. 

The following theorem of Lagrange states the order 
relations in the sequence {Cn} |22l |. 

Theorem 1 Let rational m/n be different from either 
Pi/Qi or Pi+i/Qi+i with I <n< Qi+i. Then, 

\nk -m\> \kQi ~ Pi\> \kQi+i - Pi+i \ , 

i.e., 

for any integer j satisfying 1 < j < Qi+i and j ^ Qi. 



The proof is given in Appendix B. 

Since the residence times T!„ are determined by the pa- 
rameters e„, the next corollary follows for the residence 
times. 

Corollary 1 Let 1 < j < Qi+i and j ^ Qi. 



Proof. The order relations in the sequence {£«} are 
obtained from Theorem 1. The residence time T„ is a 
monotonically decreasing function of £„, since < 0, 
as proven in Appendix C. Hence, the inequalities in 
residence times {r„} hold. □ 

Corollary 1 states that the renewal of the largest res- 
idence time r^^^ occurs when the displacement x{t) 
passes through cell Iq- {i > 1), which includes the stag- 
nation points X ~ Qi. Therefore, the largest residence 
time is given simply by 

rrnmaX rrn ^ (^i+lQi+Qi-1 

^n<Q,+ l=^Q^- ■ (l-^j 

Hence, if we know the convergent series of parameter 
fc, the positions of cells where the renewal occurs are 
completely determined. For example, for k — Lp~^, the 
stagnation points generating the largest residence time 
are determined by Eq. (jlip with qq = and qi ~ 1 {i > 

1), and given by Qi = - — — ^ — i-e., the Fibonacci 
sequence (see Line (a) in Fig. 2). 



VI. ASYMPTOTIC BEHAVIOR OF THE 
DISPLACEMENT 

According to the renewal of the largest residence time, 
the finite-time average of the velocity v{x,t) = x{t)/t is 
expected to decrease gradually as time t grows, although 
it stays positive. In this section, the asymptotic behavior 
of the displacement x{t) is investigated in detail. 



A. Asymptotic estimate of the ratio t/x{t) 

The fluctuation of the ratio t/x{t) characterizes the 
deviation from the linear increase in the displacement. 
Let us define the ratio by r{x) — t/x{t) as a function of 
the variable x under the fixed initial condition Xq = 1/2. 

It can be proven that the ratio r{x) satisfies the fol- 
lowing inequalities for large integers Qi'. 




for n^ Qi, Qi + l, ■ ■ ■ , Q^+i - 1. 



The proof is given in Appendix D by assuming the er- 
godic property of the sequence {e„}. Each term in 
Eq. (Hi]) is generally given by r (n -I- i) ~ i YTj=i 
for large integers n, which is the arithmetical average of 
the residence times up to the n-th cell. Assume that the 
dynamics of the sequence {£j}^7^ approximately er- 
godic over the interval [£Qi_i, e™^''] for large Qi where 
gmax (igjiotes the maximum of (j = 1, • • • , Q^ — 1), and 
the minimum becomes £Qi_^ from Theorem 1. Then, the 
arithmetical average can be replaced with the average cal- 
culated by the density function P{T), whose asymptotic 
form is given by Eq. PT7| . Therefore, the ratio r (Qi -I- i) 
in the left part of Eq. (|14p is approximated for large Qi 
as follows: 

~ / TP{T)dT + 

= ^lnrQ,_, +Co(fc) + ^, 

TT Qi 

where T™™ denotes the minimum residence time cor- 
responding to e"^^^, and the term Co(fc) is defined by 

Co(fc) = J^2L:' rP(T)dT-f In rQ,_,, the value of which 
depends on the behavior of P{T) in the non-asymptotic 
region. It is numerically confirmed that Co{k) converges 
to a constant for large Qi (in particular Co{(p~^) « 
0.558). From Eq. ((T3]), the ratio r (Q^ + i) can be writ- 
ten solely with information about the rational approxi- 
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mation of parameter k, 
1\ 2 



ln(Q, + Q,_i/6i,+i) 



Oi+lQi + Qi-l 



+ Ci{k), (15) 



where Ci(fc) = Co{k) — ^Imr. In the same manner, the 



ratio r [Qi+i — ^) in the right hand side of Eq. (fTi]) is 
given for large Qi+i by 



-ln(Q,+i + Q,M+2) + C\ik). (16) 



Using Eqs. (|T4|) . (|T5|l . and (|T6|) . we can estimate the ratio 
r(a;) as follows: 



I lnQ,+i + Ci(fc) < r (n + i) < I InQ, + | In (l + J^) + ^^±4i/2i + Ci(fc) 
for n = Qi, Qi + 1, • ■ • , Qj+i - 1. 

When X = n + i (n = Qi, Qi + 1, • ■ • , Qi+i ~ 1), x satisfies Qi < x < Qi+i, and Eq. (fT7|) reduces to 

I Inx + Ci(fc) < ^ < I Ina; + %i + + dik), 

for X e S",; = |n + i n = Qi, Qi + 1, . . . Qi+i — l| and its index i, 

where l/q^ < 1 and l/(gi6'i+i) < 1 are used. 



(17) 



The ratio t/x can be expressed as the sum of a logarith- 
mic term of x and some bounded function f{x), 



Ci(fc)</(x)<^ 



i ~alnx + /(x), 

2 In 2+1 



for X & Si — ^1 



n = Qi, Qi + 1, 



(18) 

Ciik) (19) 

...Q.+i-i} 



and its index i, where a = 2/7r. The above estimates 
for the ratio r{x) from (fT5)) to (fTQ]) are in good agree- 
ment with the numerical results. Figure 4 shows the 
ratios r(x) vs. x for A; = tp~^ and fc = 1/e (e = natural 
logarithm), and the lower bound for k = ip^^ given by 
alnx + Ci{(p~^). Lower bounds for each ratio are nearly 
identical, since Ci{(f~^) ~ Ci(l/e). On the other hand, 
upper bounds vary depending on the parameter k. 

If we know the value of the (i+l)-th complete quo- 
tient 9i+i and the fraction Qi/Qi-i, the lower and upper 
bounds for the ratio r{x) — t/x for x € Si are estimated 
more accurately than by Eqs. (fTS]) and ([T9|). For k = ip~^, 
relations 9i+i = (p and Qi/Qi^i c± tp are available for 
large i. Letting rL{x) and ru{x) denote the lower and 
upper bounds, respectively, the ratio r(x) is estimated as 
follows: 



'''l{x) < r(x) < ru{x) for k 
TL [x) = alnx + /l, 
ru {x) — alnx + fh + {p + p^^)/!^ 



(20) 
(21) 



where /l = a In (^1 + + Cx{^^^)- Note that ri(x) 

and ru(x) are independent of the index i. 

The inset of Fig. 4 shows that the ratio r(x) for 
k = (p~^ has zig-zag structures, which are similar to each 




X 



10000 



FIG. 4; Ratios r{x) vs. displacement x for k = ip~^ and 
= 1/e (e = natural logarithm). The dashed line is the 
common lower bound for both ratios. The inset shows the 
magnified figure of the ratio r{x) for k = (p~^ ■ The lower and 
upper bounds are given by Eqs. (|2U|I and (I2ip . respectively. 



other and appear at equal intervals in the logarithmic 
scale of X. The large peaks correspond to the stagna- 
tion points generating the largest residence time, which 



are distributed as the Fibonacci sequence, Qi 
large i. 



V5 



for 
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B. Asymptotic temporal behavior of the 
displacement x{t) 

The estimations given by Eqs. ([T8|) and indicate 
the fact that the ratio t/x{t) has an intricate time depen- 
dence. Using Eq. ([T8)) recursively, the displacement x{t) 
can be expressed as 



t 



Jn- 



a In ■ 



I In ■ 



a In ■ 



fix) 



fix) 



Provided that a In a; ^ fix), i.e., alnQ^ ^ + 
^'"^+^ + Ci(fc), is satisfied, Eq. (HHD is given by the 
t ~ ax In a;, which can be solved. Considering that the 
inverse function of y = We^ , the Lambert W function, 
is written as W{y) = -logg(- • ■log(e-„) log(e-H)(l/e)) for 
> 1 and y > e ^231, we have 



x{t) 



1 



1 



rioo 



.(i/e)' 



(22) 



where the operater S-yiu) is defined by S^yiu) 

l0g(e-«) 



ilni, 

y 



The form of the nested logarithm is not trivial. For 
k = (p~^, Eq. ([20|) . which represents the lower bound of 
t/x, can be solved as x{t) — l/[/3£^/^(l/e)] with /? = 



Here we set the n-th approximation of the lower 



bound, = l/[/3£^j/^(l/e)]. Figure 5 shows the 

ratio tjx for fc — ip^^ as a function of t with the first five 
approximations t/x''"-\t) (n = 1, 2, ■ • • ,5) : 



t 
t 



aln- + /l, 
a In - 



a In — + Jl 
a 



+ fL 



a In - 



dn- 



a In — 

a 



fL 



and so on. A higher-order logarithmic correction brings 
about better agreement with the numerical results. The 
approximation functions {a;'^"-' (t)} converge by oscillating 
to a unique function x^°°'>{t) for large t, as the order n of 
logarithmic correction increases. 



VII. SUMMARY AND DISCUSSION 

We have studied nonchaotic stagnant motion in the 
MQPGS. It was shown that the density function of resi- 
dence times obeys the inverse-square law, independent of 



+ fix) ^ 



1 i 1 1 1 1 1 1 1 

t/x 


1 1 1 1 1 1 1 1 


Approx.1 




_ Approx.2 




Approx.3 




Approx.4 




Approx.5 


iW/ •■'*'-''' 


- 


mi y ■■' 






/ 

. , 1 . , , 1 . , . 1 , , . 1 



10 



100 



1000 



10000 



t 



FIG. 5: Ratio t/x as a function of t (solid line) for k ■ 
(p~^. The other lines are approximated curves t/x^"\t) {n 
1, 2, ■ ■ ■ ,5) for the lower bound. 



the parameter k. The spatial configuration of residence 
times, however, is sensitive to k. The renewal of the 
largest residence time occurs when the displacement x(t) 
have passed through the cells Iq., i.e., by the stagnation 
points x ^ Qi, which are denominators of the conver- 
gent sequence of parameter k. Finally, the asymptotic 
behavior was given in the form of a nested logarithm. 

It is meaningful to compare the MQPGS to inter- 
mittent chaotic systems [1, H, @]- In the PM system 
mentioned in Sec. I, the residence time T of a lami- 
nar phase, starting at a reinjection point Xin scales as 
T oc cCj^ as Xin 0. Since the probability density 
function (PDF) Pmixm) of the reinjection points varies 
only slowly with Xm, i.e., Pinixm — > 0) ~ const., the 
PDF of residence times PpM iT) follows an inverse-power 
law PpAiiT) (X T"^. In particular, PpuiT) oc T^^ for 
z = 2, as in the case of Eq. pO)) . 

Confining our discussion to the Pomeau-Manneville in- 
termittency, the following two factors are essential for the 
appearance of intermittency: Local slow dynamics near 
the unstable periodic point, called intermissions, and a 
certain randomness in the seeds of residence time fluctu- 
ations. The reinjection points Xin are considered as the 
seeds in the PM system, and the parameters £„ are con- 
sidered as the seeds in MQPGS. For the two essential fac- 
tors, the local mechanisms generating intermissions are 
almost the same in the two systems. However, the seeds 
of residence time fluctuations have difi^erent properties in 
the systems. The reinjection points Xin are random in 
the sense that they are not only ergodic but also mix- 
ing (i.e., the correlation vanishes rapidly). On the other 
hand, the parameters e„ have regularity, since they are 
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only ergodic and not mixing. 

This difference is reflected in the asymptotic behav- 
ior, especially in the relaxation behavior of the long-time 
average. In the PM system, the renewal rate H{n)/n 
of chaotic bursts converges as l/lnn for z = 2, which 
is derived based on the assumption that the successive 
residence times can be considered as independent ran- 
dom variables 0, H, In the MQPGS, however, the 
occurrence rate of rapid moving phase, which coincides 
with x{t)/t, converges in the form of a nested logarithm. 
The peculiar formula reflects the quasiperiodic correla- 
tion in e„. We also note that the log-periodicity of the 
asymptotic behavior also appears in quasichaotic sys- 
tems, which exhibit weakly mixing dynamics and have 
a zero Lyapunov exponent We believe that the 

MQPGS affords another model for investigating complex 
phenomena, including slow relaxation, as well as nonsta- 
tionarity, from the viewpoint of regular system. 

In Sec. Ill we mentioned the structural instability of 
this system. The stagnation points that generate the 
largest residence time are structurally unstable, at which 
trajectories with slightly different parameter values can 
separate. This instability depends on the number- 
theoretic properties of k. We will investigate this aspect 
in a forthcoming paper [25l |. 
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Thus, 



\kQ^^P^\ 



1 



iQi + Qi-i 



Since the sequence {Qi} increases monotonically for i, 
and 9i > 1, 



\kQ,-P,\ < |fcQi-Pi| 



1 1 
< -. 



2Q1 + Qo 02qi + 1 2 



The inequality \kQi — Pi] < 1/2 assures that Pi is the 
nearest integer to kQi. As a result, eq. — \kQi — Pi\ 
holds. 



APPENDIX B: PROOF OF THEOREM 1 

This is the proof with a slight modification of the orig- 
inal one given by Lagrange |22| . Consider the equation 

kn-m = a{kQi+i - Pi+i) + (i{kQi - Pi). 

Separating the coefficient of k from the constant terms, 
we obtain two equations with two unknowns: 



n = aQj+i -I- /3Qi, m = aP^+i + (3P, 



(Bl) 



with a determinant given by Eq. (|A1[) . Since Eq. (jBip is 
transformed as (— l)*a — niQi — nPi, (— 1)'/3 = nPi+i — 
mQi+i, ct and (3 must be integers. In addition, a and /3 
cannot be zero, since m/n is different from either Pi/Qi 
or Pi^i/Qi+i. Furthermore, since n < Qi+i, a and /3 
must have opposite signs. Hence, since fcQi+i — P^+i and 
kQi — Pi also have opposite signs from Eq. (|A3p . 



\kn-m\ = \a{kQ^+i-P^+i)\ + \l3{kQ,~P^)\ > \kQ^-P,\. 
Finally, we have \kQi — Pi\ > |fcQj+i — Pi+i\ from 

Eq. ini). □ 



APPENDIX A: DERIVATION OF EQ. ([12]) 



APPENDIX C: PROOF OF THE INEQUALITY 



< 



In the theory of continued fraction, the following for- 
mulas are known: 22] 

(1) fi's and Qi's have the property that 

P,+iQ, - P,Q;+i ^ i-iy {^>-2). (Al) 

(2) The parameter k is written in another form, 

, Oi+lPi + Pi-l 



'i+lQi + Qi-1 



(*>-!)• 



(A2) 



Using Eqs. ([TT]), ((XT|) . and the difference between 

the value of k and its i-th convergent is given by 



k-^ 



%+lQi + Qi-l)Qi 



(A3) 



Let T{k, e) denote the residence time Tn{£n) defined 
by Eq. ^ abbreviating n. 



r(fc,e) 



dy 



_ 1 1 — ^ cos(27ry) — i cos(27r(£ -f ky)) 
for (fc,e) e (0,1) X (0,1/2). 



Then, the derivative of T(fc, e) with respect to e is given 
by 



dT{k,e) 



de 

where f{y) 



-47r J ^ f{y)dy, 

sin(27r(£ -|- ky)) 
(2 - cos(27ry) - cos(27r(£ + ky)))'^ 
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To know the sign of the derivative, the parameter space Thus, ^ < 0. 
is partitioned into following four regions: 

Region 1 | < e < 

f{y) > for - 1 < 2/ < 1 



Region 2 £ < min [|, iy^] 

dT{k, e) 
de 



-47r / ' f{y)dy - 47r / '^ f{y)dy, 

-4tt / fi{z)dz-4n / f2{z)dz, 
Jo Jo 

r-—— -+— 

-A-K r \h{z) + h{z))dz - An r " h{z)dz, 

Jo h-i 



where h{z) = - ^2-cos{2.ill%)-L{2.kz))- < and ^(z) = (2_eos(2^("-|l)-cos(2^fcz))^ > in ^ach domain of integra- 
tion without endpoints. In addition, 

4sm(2^-) ^m{2nkz) sm(27r^)(2 - c:os(2^A-) - (■c,s(2^-) (■c.s(2^^)) 
Ji{Z) + j2[z) ^2 _ cos(27r(2 + f )) - cos(27rfc2))2(2 - cos(27r(^ - f )) - cos(27rfc^))2 ^ ' 



Hence, ^ < 0. 
Region 3 e > max [|, 

(iT(fc,£) 



-47r / ' - 47r / ' ^ f{y)dy, 

-47r ' f3{z)dz-ATT " U{z)d 

Jo Jo 



1 1/2-e 1 1/2-e 

= -A-K r ^ {h{z) + U{z))dz - At: r " h{z)dz, 

Jo 



where /afz) = 7 ^'"/.^Z^f'f,^ ; — > 0, and f4(z) = ^TgZ^f'f,) ; — < in each domain of 

•' ^ ^ (2-cos(2-!r(2-i^^))+cos(27rfez))2 ' '' ^ (2-cos(27r(z+i^^))+cos(27rfez))2 

integration without endpoints. In addition, 

J. / X _ 4sin(27r2;)sin(27rfcz)sin(27ri^|^)(2-cos(27rfcz) -cos(27rz)cos(27ri^^)) 
(2 - cos(27r(z - i^^)) + cos(27rfcz))2(2 - cos(27r(z + i^)) + cos(27rfc2))2 

Hence, ^ < 0. 
Region 4 < e < | 

e 1/2-e 1 

dr(A;,e) 



/-k f k f2 

f{y)dy - 47r / f{y)dy - An f{y)dy, 

l_^ 1 1/2-e l-2e 1 

477 \h{z) + f2{z))dz - A-K r " {h{z) + h{z))dz - An I " ' f{y)dy. 
Jo Jo 



The first and the second terms are negative, as proven in Regions 2 and 3, respectively. The last term is negative 
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since f{y) > Q in [i - f , ^ - ^] . Thus, ^ < 0. 

As a result, the monotonically decreasing property is 
proven. □ 



APPENDIX D: PROOF OF EQ. (flil) 

For a large Qi, Eq. (fH)) is equivalent to the following 
inequalities: 



Q, -I « Qi+i-i 

, Qi+i — 1), 



(n = Qi, + 1, 



Proof of the first inequality 

The first inequality is equivalent to the following inequal- 
ity: 



Q 



> 



^^^^ 

{n = Q, + 1, 



1 



Q»+i - 1). 



(D2) 



Assume that the dynamics of {£j}"=Q■-^-l is approx- 
imately ergodic over the interval [e™™, e™**^], where 
£q.) and e"'^^^ are the minimum and maximum of 
£j = Qi + Ij ■ ■ ■ I respectively. Then, the summa- 
tion of the right-hand side of Eq. (|D2[) can be replaced 
with the average calculated by the density function P{T) 
defined by Eq. ^TU\\ as follows: 



1 



n- Qi 



j = Q^ + l 



TP{T)dT < / TP(T)dT, 
(D3) 



where the last inequality comes from T™'^'^ < Tq., due 
to Corollary 1. The last integral in Eq. (|D3p is given by 



TP{T)dT ^ - \n{Q^+Q^-l/e^+^)+- ln0,+i+Ci(fc). 

(D4) 

Using Eq. (|D4)) and Eq. ((151) for the left-right hand side 
of Eq. jnil), we find 



1 



Qi 



. . (e.+i -21n6',+i)Q, -hQi_i 



E^. 



TrQi 



where 6* — 2 In > is used. Therefore, Eq. (|D2|) . i.e. 
the first inequality holds. 

Proof of the second inequality 

The middle term in Eq. (|Dip can be written as 



n 



Qi+i-i Q,+i-i 

E E 



_ i=l j=n-|-l 

Since the summation X]j=in+i ^-Iso does not have any 
term larger than Tg^, due to Corollary 1, it is estimated 
as 



I 



where the ergodicity of {ejjjll^Yi assumed again. 
Hence, 



1^ ^ 1 



(Q,:+i - 1) / TP{T)dT - (Q,+i -1-72)1"' TP{T)dT 



/ TP{T)dT^- 



Thus, the second inequality holds. 



□ 



[1] P. Berge, M. Dubois, P. Manneville, and Y. Pomeau, J. [3] W. J. Yeh and Y. H. Kao, Appl. Phys. Lett. 42, 

Phys. (Paris), Lett. 41, L341 (1980). 299 (1983); H. Seifert, Phys. Lett. 98A, 213 (1983); 

[2] Y. Pomeau, J. C. Roux, A. Rossi, S. Bachelart, and C. Vi- Ch. Noldeke and H. Seifert, Phys. Lett. 109A, 401 (1985) 

dal, J. Phys. (Paris), Lett. 42, L271 (1981). [4] Y. Pomeau and P. Manneville, Commun. Math. Phys. 



74, 189 (1980); P. Manneville, J. Phys. (Paris), Lett. 41, 

1235 (1980). 

[5] T. Gciscl and S. Thomae, Phys. Rev. Lett. 52, 1936 
(1984). 

[6] P. Gaspard and X.-J. Wang, Proc. Natl. Acad. Sci. USA 

85, 4591 (1988). 
[7] Y. Aizawa, Chaos, Soliton and Fractals 11, 263 (2000). 
[8] T. Akimoto and Y. Aizawa, Nonlinear Phenomena in 

Complex Systems 9, 178 (2006). 
[9] D. R. Cox, Renewal Theory (Mcthuen, New York, 1962). 
[10] A. Prasad, V. Mehra, and R. Ramaswamy, Phys. Rev. 

Lett. 79, 4127 (1997). 
[11] A. Witt, U. Feudel, and A. Pikovsky, Physica D 109, 180 

(1997). 

[12] A. Venkatesan, K. Murali, and M. Lakshmanan, Phys. 

Lett. A 259, 246 (1999). 
[13] A. Venkatesan, M. Lakshmanan, A. Prasad, and R. Rar 

maswamy, Phys. Rev. E 61, 3641 (2000). 
[14] S.-Y. Kim, W. Lim, and E. Ott, Phys. Rev. E 67, 056203 

(2003). 

[15] S. H. Strogatz, Nonlinear Dynamics and Chaos: with Ap- 



11 

plications in Physics, Biology, Chemistry, and Engineer- 
ing (Addison- Wesley, Reading, MA, 1994). 

[16] A. T. Winfroo, The Geometry of Biological Time, 2nd ed. 
(Springer- Verlag, Berlin, 2001). 

[17] I. Zapata, R. Bartussek, F. Sols, and P. Hanggi, Phys. 
Rev. Lett. 77, 2292 (1996). 

[18] W.C. Stewart, Appl. Phys. Lett. 12, 277 (1968). 

[19] D.E. McCumbcr, J. Appl. Phys. 39, 3113 (1968). 

[20] A. Barone and G. Paterno, Physics and Applications of 
the Josephson Effect (John Wiley & Sons, New York, 
1982). 

[21] V. I. Arnold and A. Avez , Ergodic problems of Classical 

Mechanics (Benjamin, New York, 1968). 
[22] A. M. Rockett and P. Szusz, Continued Fractions (World 

Scientific, Now York, 1992). 
[23] Y. Clio and K. Park, The American Mathematical 

Monthly, 108, 963 (2001). 
[24] G. M. Zaslavsky. Phys. Rep. 371, 461 (2002). 
[25] T. Mitsui, Prog. Theor. Phys. Suppl. 173, 243 (2008). 



